Constraining Cosmic Expansion and Gravity with Galaxy Redshift Surveys 
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We analyze the science reach of a next generation galaxy redshift survey such as BigBOSS to 
fit simultaneously for time varying dark energy equation of state and time- and scale-dependent 
gravity. The simultaneous fit avoids potential bias from assuming ACDM expansion or general 
relativity and leads to only modest degradation in constraints. Galaxy bias, fit freely in redshift 
bins, is self calibrated by spectroscopic measurements of redshift space distortions and causes little 
impact. The combination of galaxy redshift, cosmic microwave background, and supernova distance 
data can deliver 5-10% constraints on 6 model independent modified gravity quantities. 



I. INTRODUCTION 

The power spectrum of matter density perturbations 
contains information on both the cosmic expansion his- 
tory and growth history. Comparison of these fundamen- 
tal evolutions not only tightens constraints on the cosmo- 
logical model but enables tests of the framework, such as 
the validity of general relativity as the theory of gravity. 
However, most cosmological parameter estimation, from 
cither current or projected data, either assumes general 
relativity, and hence fixes the growth history to be de- 
termined by the expansion history, or assumes a ACDM 
expansion history to test gravity. 

If either assumption is incorrect then the result will be 
biased. Thus, even if one is only interested in constrain- 
ing the expansion history and dark energy properties, for 
example, one still must fit for gravity in order to obtain 
robust results for the expansion. Current data cannot 
tightly constrain either quantity without ad hoc assump- 
tions such as no time dependence. The volume and red- 
shift reach of cosmological surveys is rapidly improving, 
however, and projections for the next generation are for 
the ~ 10% accuracy level on the time variation of dark 
energy properties or gravity. 

The question we address is how the constraints are 
affected when the matter power spectrum data mea- 
sured by galaxy redshift surveys are fit simultaneously 
for varying dark energy equation of state, scale- and 
time-dependent gravity, and astrophysical evolution such 
as galaxy bias. Aspects of this have been consid- 
ered in the literature with some restrictions: for exam- 
ple [1] parametrized gravity with the scale- and time- 
independent gravitational growth index 7 and a single 
bias parameter, while [2] extended this to binned bias; 
[3, 4] assumed scale-independent gravity with a particu- 
lar time variation but allows scale- and time-dependent 
bias; [5] applied a motivated scale- and time-dependent 
gravity parametrization to current data, restricting to 
quantities independent of galaxy bias. (Note that more 
general simultaneous fitting, such as with principal com- 
ponent analysis, has been applied to other cosmological 
probes such as weak lensing, e.g. [6, 7].) 



With the aim of deriving more general results that 
avoid assuming a particular model, for expansion his- 
tory, gravity, or astrophysical bias, we use scale- and 
time-dependent bins of gravity and time varying effec- 
tive dark energy equation of state and bins of galaxy 
bias. Specifically, we fit for gravitational modifications 
to the two Poisson equations (for the behavior of mat- 
ter and of light) as freely floating values at high and low 
redshift and large and small scales, i.e. in bins of redshift 
z and wavenumber k. Galaxy bias is fit as independent 
values in 17 bins of redshift. Dark energy evolution is 
treated through the highly accurate WQ-w a parametriza- 
tion. Due to the presence of non-Gaussian covariances 
between the many parameters we carry out a Markov 
Chain Monte Carlo exploration of parameter space using 
simulated next generation cosmological survey data. 

In Sec. II we explain the treatment of dark energy, 
gravity, and galaxy bias in detail. The simulated data 
sets used are presented in Sec. III. We analyze the results, 
with particular attention to covariances and the effect on 
constraints of not assuming fixed expansion, fixed grav- 
ity, or fixed galaxy bias in Sec. IV. Implications for the 
science reach of next generation surveys are discussed in 
Sec. V. 



II. DARK ENERGY, GRAVITY, AND 
GALAXIES 



Matter density perturbations grow under gravitational 
instability but at a rate suppressed by cosmic expansion. 
Therefore the evolution of large scale structure clustering 
depends on the matter density fi m , the expansion rate 
H, and the laws of gravity. In addition, galaxy redshift 
surveys do not directly measure the mass power spec- 
trum in real space but the galaxy clustering in redshift 
space. This introduces two additional ingredients: the 
statistical distribution of galaxies may be biased relative 
to the mass (i.e. dark matter), and the velocity field of 
the galaxies, generated by the gravitational potentials of 
the structures, adds anisotropic distortions to statisti- 
cally isotropic density field. 
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The rcdshift space distortion depends on the angle with 
respect to the line of sight and its overall amplitude is 
determined by the growth rate of structure at the given 
redshift, f(z). This carries with it important additional 
cosmological (and gravitational) information, enhancing 
the reach of spectroscopic galaxy surveys. Galaxy bias 
6, however, has the potential to confuse extraction of the 
total amplitude of growth up to the given redshift. Since 
it enters without angular dependence, one can fit sepa- 
rately the two effects, although each is convolved with 
the overall mass fluctuation amplitude often denoted by 
eg (z). That is, cosmological information comes in the 
form of bas(z) and /crg(z) [8]. The amplitude ag,(z) is 
proportional to the linear growth factor D(z), which is 
related to the gravitational potential through a Poisson 
equation, and so this growth history depends on both the 
expansion history and gravity. If additional data directly 
related to the mass fluctuations (rather than galaxies) is 
available, for example through weak gravitational lens- 
ing, then D(z) can be separately determined, although 
lensing also depends on gravity, in its own way. 

Thus to draw general robust conclusions on expan- 
sion and dark energy properties, for example, we should 
fit not only for H(z) but simultaneously for gravity - 
the two histories then determining D(z) and f(z) = 
— d\nD/d\n(l + z) - and galaxy bias. Moreover, the 
way we treat each of these quantities should be valid 
over a wide range of cosmologies so that the results are 
not overly model dependent. 

For the expansion history we work within the flat 
Friedmann-Robertson- Walker framework, where the ex- 
pansion rate is described by the matter density f2 rn as 
a fraction of the critical density and the time dependent 
dark energy equation of state w(z) = wq + w a z/(l + z). 
This form for w(z), obtained from exact solutions for a 
wide range of dark energy and gravity models [9], has 
been demonstrated as a calibration relation with an ac- 
curacy of 0.1% for obscrvables [10]. Perturbations in the 
effective dark energy (i.e. even if there is no physical dark 
energy, just a modification of gravity) are treated consis- 
tently within the equations of motion. 

Note that while many papers take an exact ACDM ex- 
pansion history when allowing for modified gravity, this 
is highly model dependent. While simple f(R) gravity 
models can be treated in terms of the scalaron mass, 
which determines both the growth and expansion be- 
haviors so that effective screening within the solar sys- 
tem (really, small df/dR) requires negligible deviations 
of w(z) from —1, this is not generally true. For example, 
DGP gravity [11] has strong deviations in both growth 
and expansion from ACDM, though as it is a one pa- 
rameter model they are tied together. Galileon models, 
however, have sufficient freedom that the significant de- 
viations in expansion history can appear alongside mod- 
erate growth deviations (see, e.g., [12]). Thus fitting for 
expansion and gravity independently is most generic and 
free from possibly unwarranted assumptions. 

To account for possible extensions beyond general rel- 



ativity, and their effects on matter growth and lensing of 
light, we solve for the two metric potentials <p an d 4> m 
modified Poisson equations 

V 2 V> = inG N a 2 8p x G ma tter (1) 
V 2 (0 + r/>) = 8nG N a 2 5p x G Iight . (2) 

The rigor and completeness of these equations, together 
with those for the matter density and velocity from 
energy-momentum conservation, is discussed in detail in, 
e.g., [13-16] (though here we have simplified the nota- 
tion) and clearly related to the nonrelativistic and rela- 
tivistic geodesic equations in [17]. 

Paramctrization of the generally time- and scale- 
dependent functions G ma tt e r(&; z) and Gii g ht(fc, z) can be 
specialized to forms expected from scalar-tensor or DGP 
gravity, for example (see [5] for a unification of the two), 
or kept free form through a principal component analysis 
[6, 7, 15, 18], say. Since we wish to carry out a substan- 
tially model independent analysis, but be able to inter- 
pret clearly physically the constraints on gravity, we fol- 
low the high vs. low redshift, large vs. small scale binning 
approach of [13, 16]. This allows for 8 free gravity param- 
eters: Gmatter and Giight values in bins from z = 0-1 and 
z = 1-2 (for z > 2 their values are set to 1, i.e. behaving 
as general relativity in the high redshift, high curvature 
universe), and in wavenumber bins 10~ 4 Mpc _ < k < 
10- 2 Mpc _1 and 0.01 Mpc" 1 < k < 0.1 Mpc" 1 . Because 
of the uncertainty in the appropriate growth equations in 
the nonlinear density region at higher wavenumbers for 
an arbitrary gravity theory (e.g. the presence of various 
screening mechanisms [19]), and in the velocity induced 
redshift space distortions to the galaxy power spectrum, 
we conservatively do not use data at k > 0.1 Mpc -1 . 

Galaxy bias is treated through 17 free parameters rep- 
resenting the values b(z) in bins of width 0.1 in redshift 
from z = 0.1-1.8. Since we only include large scale 
information we take the bias to be scale independent. 
We phrase the bias amplitude in terms of the comov- 
ing clustering expectation, b(z) = bo(z) D(z = 0)/D(z) 
and fit for bo in each redshift bin. We use independent 
sets {&LRG( z i)} an d {^elg( z 0} f° r bo(z) for luminous red 
galaxies and emission line galaxies in the survey. 

In addition to these parameters we also fit for the 
physical baryon density Qbh 2 , physical matter density 
il m /i 2 , ratio of cosmic microwave background (CMB) 
sound horizon to angular diameter distance to last scat- 
tering 9, optical depth to reionization r, scalar spectral 
tilt n s , amplitude of primordial scalar perturbations A s , 
and any astrophysical nuisance parameters appropriate 
to the data sets. 



III. GALAXY REDSHIFT SURVEY SCIENCE 

Spectroscopic galaxy surveys provide three dimen- 
sional information on galaxy clustering and direct mea- 
surements of the galaxy density and velocity fields, thus 
probing the growth history. BigBOSS [20] is planned to 
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survey 50 (Gpc//i) 3 of cosmic volume, obtaining redshifts 
of 20 million galaxies from z ~ 0.5-1.8. (In addition, it 
provides data on quasar clustering and on neutral hydro- 
gen density fluctuations through the Lyman alpha forest, 
ranging over z = 1.8-3.5, but conservatively we do not 
include constraints from those probes.) 

The redshift space galaxy power spectrum to be com- 
pared with observations is related to the linear theory 
mass power spectrum calculated from the equations of 
motion (Poisson equations and energy-momentum equa- 
tions) by 

p*(k, z ) = (b+fv 2 ) 2 p,; nMn (k,z) , (3) 

where the squared factor is the Kaiser correction [21]. 
This is a good approximation strictly in the linear regime; 
since we use only k < 0.1 Mpc -1 this is not too unreason- 
able. How to map the isotropic mass power spectrum to 
the anisotropic galaxy power spectrum is not settled be- 
yond linear theory so we do not attempt to go beyond the 
Kaiser formula here. (See [2, 22] for explorations of be- 
yond linear theory constraints on the scale independent 
gravity parameter 7 as a function of maximum wavenum- 
ber, and [23] for specific f(R) gravity models.) 

Uncertainties on the measurement of the galaxy power 
spectrum depend on the volume surveyed in each redshift 
shell and the number density n of each type of galaxies 
used to sample the density field. The volume is deter- 
mined by the area of sky (14000 deg 2 for BigBOSS) and 
the cosmological distance-redshift relation. We use the 
redshift distributions n(z) for the two types of galaxies 
given in [20]. See [2, 24, 25] for further details on calcu- 
lating the likelihood with galaxy power spectrum mea- 
surements. 

In addition, we include simulated Planck CMB data. 
This helps to constrain the cosmological parameters, es- 
pecially the primordial ones, and the gravitational mod- 
ifications through the integrated Sachs- Wolfe effect and 
CMB lensing. The ISW effect has sensitivity to Gn g ht 
since this Poisson equation governs the behavior of rel- 
ativistic geodesies; CMB lensing is sensitive to both 
G ma tter through growth in the matter power spectrum 
and Gnght through light deflection. We include CMB 
lensing through the deflection angle power spectrum as 
prescribed by [26]. (See [27] for its use in constraints 
on the gravity parameter 7.) For constraints on the ex- 
pansion history we employ simulated supernova distance- 
redshift measurements of 1800 supernovae over z = 0-1.5 
with a systematic floor of 0.02(l+z)/2.5 magnitudes, and 
include the supernova absolute magnitude parameter A4 
as a nuisance parameter. 

IV. CONSTRAINING EXPANSION AND 
GRAVITY SIMULTANEOUSLY 

We carry out a Markov Chain Monte Carlo analysis 
to find the joint posterior likelihoods of the parameters, 
with our modified versions of CAMB [28] and CosmoMC 



[29]. Our baseline analysis consists of projected galaxy 
clustering, CMB, and supernova data, allowing for inde- 
pendent fits of the 8 binned gravity parameters Gi; g ht and 
Gmattcr, 27 binned galaxy bias parameters for LRG and 
ELG, the dark energy equation of state parameters Wq, 
w a , plus the other cosmological and nuisance parameters. 
To speed the convergence of our Markov Chain Monte 
Carlo implementation, we fit for the maximum-likelihood 
bias combination at each step in cosmological parameter 
space. Because of this, our quoted credibility regions 
should be considered conservative, as outlying points in 
cosmological parameter space are assigned a larger pos- 
terior density than they might be if we marginalized over 
the bias parameters. We then consider the following vari- 
ations to explore the impacts of physical ingredients: 

• Expansion: fix to ACDM background, i.e. wq = — 1, 
w a =0 

• Gravity: fix to general relativity, i.e. G ma ttcr = 1, 
Giight = 1 for all bins 

• Galaxy bias: fix to inverse growth, but fit for con- 
stant b ELG 

We emphasize that our baseline case varies all these in- 
gredients. 

The results for the gravitational sector in the base- 
line case and its physical variations are shown in Figs. 1. 
We see that for next generation data sets there is little 
degradation in fitting for dark energy equation of state 
and galaxy bias simultaneously with gravity, giving great 
promise to our ability to test the cosmological framework. 
Adding parameters for non-ACDM expansion has only a 
percent level effect (except ~ 15% on G ma tt CT at large k 
and z) , while bias parameters only affect G ma ttor at large 
k and only loosen the constraints by at most ~ 50%. 

The quantity Gn g ht will be determined to ~ 5% (at 
68% CL) independently at low and high redshift, and 
small and large wavenumber, while G ma tter will be tested 
to ~ 8% on small scales (k > O.OlMpc" 1 ) where the 
galaxy survey data is most important, and to ~ 35% on 
large scales (k < O.OlMpc -1 ). The much weaker con- 
straints at large scales, where we assume no galaxy data, 
demonstrate the importance of the galaxy power spec- 
trum measurement in constraining gravity. Even so, the 
constraints on large scale G ma tter show more than a fac- 
tor of two improvement over the yellow shaded contours 
in Figure 11 of [16] (their V is our G ma tter), due to our 
inclusion of CMB lensing. 

The gravity parameters have little covariance with 
each other (including between large and small scale 
bins of the same parameter, not shown), except mildly 
at low redshift. Excellent complementarity exists be- 
tween the quantities and the probes: galaxy data con- 
strains G ma ttcr, CMB data predominantly constrains 
Gnght (again see Figure 11 of [16] and Figure 7 of [13]), 
and supernova data constrains the expansion dynamics 
of wq, w a , ignoring the gravity parameters. 
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FIG. 1. 68% and 95% projected confidence limit contours for the gravitational parameters binned in low and high redshift, low 
and high k are plotted for the baseline case of fitting for gravity, expansion, and galaxy bias (black curves), and the variations 
of fixing expansion by w = — 1 (light green curves) or fixing galaxy bias to inverse growth (dark red curves). 



The expansion sector is also not strongly affected by 
the different cases, when next generation data is avail- 
able. Figure 2 shows the contours in the baseline case, 
when fixing galaxy bias, and when fixing to general rel- 
ativity. The WQ-w a contour increases in area by ~ 30% 
when including the fit for the 8 gravity parameters, with 
~ 2% loosening in wo and ~ 20% in w a determination. 
As expected from these and the previous results, plots 
of the crosscorrelations of gravity and expansion param- 
eters (not shown) display little covariance, with only a 
minor correlation between w a and G ma ttcr at large k. 



V. CONCLUSIONS 

Testing gravity on cosmic scales is crucial to under- 
standing acceleration and fundamental laws of physics. 
The inability of present data sets to significantly con- 
strain the laws of gravity affecting matter (i.e. G ma tter) 



is the largest obstacle to efforts to test general relativity 
on cosmological scales and to constrain modified grav- 
ity theories. Next generation astronomical surveys pro- 
pose to remove this obstacle by presenting us with de- 
tailed redshift-spacc maps of the distributions of galaxies 
throughout the Universe. 

We have shown that these maps will, indeed, be able to 
deliver on their promise to constrain gravity at cosmolog- 
ical scales, even given our imperfect knowledge of both 
the cosmic expansion history and the bias factor relat- 
ing galaxy to dark matter clustering. While a theoretical 
framework predicting the galaxy-dark matter bias would 
certainly be welcome, the results show it is not requisite 
for cosmological tests of gravity. Physics such as a gener- 
alized expansion history and astrophysics such as galaxy 
bias should be included in the fits (rather than assumed 
known) to avoid erroneous conclusions about gravity, and 
can be included without fear of degrading observational 
constraints beyond usefulness. Conversely, future exper- 
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used here avoids restriction to a specific model, and 
has excellent orthogonality among its quantities and be- 
tween its quantities and the expansion and galaxy bias 
parameters. Strong complementarity exists as well be- 
tween the cosmic probes: galaxy surveys testing G ma tter, 
CMB measurements constraining Gu g ht, and supernova 
distances measuring expansion. Together, the next gen- 
eration combination of these data will deliver 5-10% mea- 
surements of 6 gravity quantities, plus 35% measures of 
the remaining 2. These can potentially be improved fur- 
ther by going beyond the purely linear regime (we con- 
servatively cut off k > 0.1 Mpc -1 ) once modified gravity 
effects on the redshift-space galaxy power spectrum are 
better understood, by higher quality CMB lensing from 
ground based polarization experiments, or by inclusion 
of crosscorrelations between galaxy maps and the CMB. 



FIG. 2. 68% and 95% confidence limit contours for the effec- 
tive dark energy equation of state parameters wo and w a are 
plotted for the baseline case of fitting for gravity, expansion, 
and galaxy bias (black curves), and the variations of fixing to 
general relativity (light blue curves) or fixing galaxy bias to 
inverse growth (dark red curves). 

iments can proceed without fear of their science output 
being obscured by our ignorance. 

The model independent approach to testing gravity 
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